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It is shown that the entanglement-structure of 3- and 4-qubit states can be characterized by 
optimized operators of the Mermin-Klyshko type. It is possible to discriminate between pure 2- 
qubit entanglements and higher entanglements. A comparison with a global entanglement measure 
and the i-concurrence is made. 
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I. INTRODUCTION 



There seems to be no doubt in the literature that en- 
tanglement of quantum mechanical states is one of the 
important ingredients in the broad field of quantum in- 
formation theory. Many protocols in this area are based 
on entangled states, it is the basis of quantum cryptog- 
raphy 1], super dense coding 2], teleportation ^ and 
other fields. Although it is an important ingredient en- 
tanglement is even called puzzling, and especially genuine 
multipartite entanglement for 3- and 4-qubits is a field 
of active research, because only for two qubits a correct 
measure for their entanglement is available. 
In this paper we will give a new method of characteriza- 
tion of 3- and 4-qubit entanglement. 
In the following we will show how different measures for 
quantifying entanglement can be applied to a model spin 
system which can be used as a basis for many differ- 
ent experimental setups. It is a one-dimensional Heisen- 
berg spin system with different arrangements for 3- and 
4 spins. We will compare a global entanglement measure 
U with the results of Bell inequalities in the form pro- 
posed by Mermin and Klyshko 5, 6, 7] which means that 
we look for a measure with optimized polynomial spin op- 
erators. It will be shown that the optimized polynomials 
measure the different degrees of entanglement astonish- 
ing well. 

The paper is structured as follows: In chapter II a special 
form of Bell inequalities is described for 3- and 4-qubit 
systems. We also define a new method how to handle these 
polynomial inequalities. In chapter III we shortly dis- 
cuss the global entanglement measure and describe its 
relation to concurrences as well as to 3-qubit tangle mea- 
sure. In chapter IV these measures are applied to the 
general Greenberger Home Zeilinger (GHZ) state for 3- 
qubits and the surprising result is that the optimized 
Mermin-Klyshko operators are well suited to describe the 
entanglement as well as the global entanglement over a 
wide range of parameters. In chapter V we now apply 
these different measures to the special form of a 3-qubit 



Heisenberg spin system and again it turned out that for 
the pure eigenstates and for the superposition of these 
eigenstates the optimized Bell operator as well as the 
global entanglement measure are in excellent agreement 
describing the entanglement. In chapter VI we have a 
look at a 4-qubit system and compare it with the 3-qubit 
results. A discussion follows in the concluding chapter 
VII. 



II. SPIN-POLYNOMIALS FOR 

ENTANGLEMENT MEASURE 

From Bell-Inequalities various polynomials are known 
for entanglement-classifications. We use here the 
Mermin-Klyshko polynomials which were discussed by 
Yu et al. 8] in their classififcation scheme. We propose 
a special optimization procedure to give a quantitative 
analyse of Heisenberg spin models. 

We test our method for 3-qubit systems and then ap- 
ply it to 4-qubits, where it is shown that the resulting 
optimized Mermin-Klyshko operators are able to quanti- 
tatively describe the total entanglement measures. Fur- 
thermore, comparing it with the sum of the 2-qubit con- 
currences the difference gives a measure to the additional 
3- and 4-qubit entanglements. We note already here that 
the optimization procedure yields many equivalent min- 
ima, therefore it is not possible to extract directly a single 
well defined operator polynomial. 



A. 3 qubits 

For 3-qubits the polynomials can be written as prod- 
ucts of spin operators 8] : 

Fg = {AB' + A'B)C + {AB - A'B')C' (1) 

f;^ = {AB' + A'B)C' - [AB - A' B')C (2) 

where the operators are written as sums of Pauli matri- 



Ai') 



d^'^-<7A, 



B^'^^P^-ah, CW=c<') 



c^c 
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with o^'', b^'^ and c^'' normalised vectors and the Pauli 
matrices a a, (Jb, f?Cj refering to the qubits A, B and C, 



withCT, = (crf,0-f,0-f). 

The classification for pure 3-qubit states is as follows. 
(For details we refer to |3|-) We look for the maximum 
of the absolute values of the expectation values of these 
polynomials and find for pure product states: 



4l(^3)p|,|(i^:^)p|}<2 



(3) 



The other inequalities found in 8] can be written as 

(F^)l + (F^)l<2' (4) 

if the state is 2-qubit entangled and 

{F,)l + {Ff,)l<2^ (5) 

if the state is 3-qubit entangled. 

Our method consists in a numerical optimization of the 
components of the vectors a^'' , 6^'^ and c^'^ of these poly- 
nomials. We used the NAG library function e04ucc [13 
with randomly chosen initial conditions. 
Since it was not clear from [^ we applied two different 
methods. In the first approach we look for the expec- 
tation value of F^ and maximize it. The polynomial 
F^ is calculated with these paramters and then the sum 
(Fs)^ + (F^)^ is determined. 

In the second approach we directly maximized the sum 
of the squares since all these inequalities are sufficient 
but not necessary. By comparing the results of the op- 
timization with each other and with other measures of 
entanglement we found that the first described method, 
the optimization of the F3— operator yields the best in- 
formation. 

Since we have not seen this numerical investigation even 
for the simplest states used in the literature we cite 
the following results, obtained with the first described 
method. 
For the GHZ-state 0] (|000) + |111))/V2: 



max{|(^3)GHz|} = 4.00 
^m)GHz\=0: 



(F 



3/ GHZ 



+ (^3; 



GHZ 



= 16.00 



while for the so-called W-state [T^ (|001> + |010) + 
|100))/\/3 we can list the following results which show 
that there seems to be a 3-party entanglement although 
compared to GHZ it has not the maximum possible value. 



max{|(i^3)w|} = 3.05 




^|(F^)v^|=0.05 {F^f^ + iF^fw 


= 9.305 


B. 4 qubits 





The sufficient conditions for 4-qubits are a little more 
involved since we have to introduce an additional spin 
polynomial for the 4**^ qubit but we can write altogether 



Fi = -(D 

F' = -(D 

^ 2^ 



D')®F3 + ^{D-D')< 



F' 
^3 



D')®F;^ + -{D' -D)®F:i 



(6) 

(7) 



where F3 and F^ are defined in Q resp. (0) and 
£)(') = (?') • ah are the spin operators on the 4**^ qubit D. 
The classification scheme after Yu et al. Q is as follows. 
Product states fulfil the following inequality: 



4l(^4),U(Fi)p|}<2 



(8) 



And for distinguishing different kinds of entanglement 
one can use the following scheme which gives sufficient 
but not necessary classification. 



2-qubit entanglement: (i<4 
3-qubit entanglement: {F^^ 
4-qubit entanglement: (-F!^ 



{Fi 



^< 



{Fi)l < 16 



{Fi)l < 32, 



where the description j^ is as follows. 4-qubit entangle- 
ment means a state with fully entangled 4-qubits, 3-qubit 
entanglement describes a product-state of one qubit with 
fully entangled 3-qubits, and 2-qubit entanglement can 
be a product of two 2-qubit entangled states or a 2-qubit 
entangled state as product with two single qubits. 
Before applying these inequalities to the spin-systems we 
discuss another useful measure. 



III. I-CONCURRENCES AND GLOBAL 

ENTANGLEMENT 

The original measure Q of a many qubit pure state \ip) 
was introduced by Meyer and Wallach Q. It was later 
shown by Brennen [ij that this kind of global entangle- 
ment can be written as 



Qim^2[l-l/nJ2Tr{pl)], 



(9) 



k=l 



with pfc, the density matrix reduced to a single qubit 
k. It is interesting to note that there can be introduced 
the so-called i-concurrence [13 which also is directly re- 
lated to the reduced density matrix pA of a subsystem 
A. This i-concurrence measures the entanglement be- 
tween two subsystems A and B and can be written as 



ICa^b = V2[l 



Tr{p\)] 



(10) 



In the following we use the notation ICa-b = ICa- 
We find as first result that the global entanglement Q 
is directly related to the sum of the squares of the i- 
concurrences of the 1-qubit subsytems of a A^ qubit state 



N 



q = i/nY,ic? 



(11) 



A. 3 qubits 

For the special case of 3-qubits one can introduce the 
so-caUed tangle T123 [13 which in a sense describes those 
contributions to the i-concurrences which are not de- 
scribed by 2-qubit concurrences |14.[l5l| 



IC'i = Cf^ 
IC2 — C12 



ici 



^2 



Cl3 
C23 

a 



23 



We can sum these relations up 



3 



c 



T123 
T123 
T123 



'^23) +3ti23 



(12) 
(13) 
(14) 



(15) 



and introduce this into the global entanglement. It is 
nicely seen that for 3 qubits this consists of the sum of 
squared 2-qubit concurrences plus the additional tangle: 



Q^^\ 



^13 



^23) 



'''123 



(16) 



The total entanglement measure is the sum of different 
entanglement contributions. 



B. 4 qubits 

These nice results for 3-qubits cannot easily be ex- 
tended to 4-qubits since there is no equivalent definition 
of the corresponding higher tangle. But to give an im- 
pression of the power of the description with a global 
measure one can look for special qubit states were there 
are effectively only 2-qubit concurrences. 

|(/)) = aillOOO) -I- aalOlOO) -f a3|0010) -I- a4|0001) (17) 

One easily finds that the i-concurrences are sums of 2- 
qubits concurrences and therefore the global entangle- 
ment can be written as 



^W, 



Ci3 



^■14 + C23 + C24 + C34) . (18) 



Again this indicates a good total measure of entangle- 
ment by the value of Q. 



IV. APPLICATION TO GENERALIZED 
GHZ-STATE 

One result of our investigations is that the comparison 
of sufficient conditions from the Bell inequalities and the 
global expression Q is an appropriate measure for the 
entanglement of 3- and 4-qubits. 

As a first test we consider the generalised GHZ state for 
3 qubits written as 




FIG. 1: Optimization of the expectation value of F3 for 
the 3-qubit generalized GHZ-state as a function of g^\ the 
course of niax{|{_F3)|} is fitted with y = Cog-y/l — <7^, with 
Co = 8.01; the course of (Fs)^ + {F^f with y = cig^(l - g^) 
and ci = 64.01; maxK^s)^ -I- {F^)^} marks the optimization 
of the squared inequalities. The dashed lines mark the bound- 
aries arising from the inequalities resp. (QJ. See text for 
more discussion. 



with g G [0, 1]. It is well known that there are no 2-qubit 
concurrences so that the remaining i-concurrences 



ICi = IC2 = IC3 = '2g^/l~, 



(20) 



mainly measure the tangle of the state which is of course 
parameter dependent and from our formula H16() it can 
be seen that Q just measures this tangle: 



Ti23^4:g^{l~g^) 



(21) 



5|000> + Vl-.9'|111> 



(19) 



In fig. ^the results of the Bell optimization are plotted 
as a function of g^. max{|(i^3)|} means the optimiza- 
tion of the expectation value of F3. Then the resulting 
parameters are introduced in the expectation value of 
F^ (Notation: K-Fg)!) and the squared inequalities Q) re- 
spectively © (Notation: (Fg)^-!- (^3)^ ). This procedure 
yields results that are nearly identical to the parameter 
dependence of the i-concurrence and the tangle resp. Q 
(cp. (|2Hl 'l over a large parameter range. This fact is 
underlined with the fit of the results of the inequalities 
by the results of the i-concurrences resp. the tangle. We 
therfore conclude that our maximization procedure yields 
the correct information for the entanglement. (Note how- 
ever that the optimization procedure as with all search- 
algorithms may fail to find the appropriate maximum.) 
Only for a small parameter range (0 < g < 0.39 and 
0.92 < g < 1) the inequalities are not sufficient compared 
to the calculated tangle. These results show that at the 
boundaries of the parameter values the optimization of 
F3 could have problems but over a large range of the pa- 
rameter values the optimized F3 gives a good measure 
of entanglement as the tangle itself although we have no 
direct proof of the equivalence of these two meassures. 
Quite remarkable is here the fact that the optimiztion of 
{Fsf + (F^)2 (Notation: max{(F3)2 + {F^f}) yields for 
the parameter values 0.39 < g < 0.92 no sufficient crite- 
rion for entanglement. 

The points {g = 0.26 and g = 0.97) where the inequal- 
ity ||2Il is not violated agree with the results derived by 



Scarani and Gisin H^. 

This successful description of the generalized GHZ-states 
encourages us to describe the entanglement of more com- 
plex systems. 



was seen for the generalized GHZ-state the parameter de- 
pendence of the states gives insight into the efectiveness 
of different entanglement measures. 



V. PURE 3-QUBIT HEISENBERG STATES 

As in reference [l3| already discussed, Heisenberg spin 
sytems are good models for various experimental realiza- 
tions of multi-qubit systems. Here we look for a special 
chain of 3 qubits where the interaction between qubit 
1 and 2 and 2 and 3 are given by a certain interaction 
constant while between 1 and 3 we have doubled the in- 
teraction. Our main purpose is, however, to study the 
anisotropy effect of this Hamiltonian given by 



4 



.y^y 



^2"3 



ay,ay, + IS.iji,a 



2"3 



{a'la'l -f a\a\ + I\gIg\ 



(22) 



with the anisotropy coefficient A. The eigensystcm of 
this Hamiltonian is calculated in the computational basis. 
The eigenvalues and eigenstates are given in Table I. 

The eigenstates are of course partially degenerated be- 
cause of the spin symmetry of the system. This can be 
easily lifted by an applied field 



H' = H 



M 






(23) 



so that in the following we think of the different eigen- 
states as pure states and discuss only the parameter de- 
pendent eigenstates IV's) to 1-08) because we are interested 
in the change of entanglement with different anisotropy 
strengths. Since the 2-qubit measures are known we will 
get more insight into pure 3-qubit entanglement. As it 



A. The states Jt/js), l-^e) 

We start with the states jV's) and \tpe) which yield the 
same results in measuring the entanglement. The con- 
currences are calculated to be: 



Cl2 — C23 — 



C'la — 



r, r;2-(A-2)7?- 

In the limit A — > cxd the states have an easy form 
1 

1 



(24) 



m 
1^6) 



:(|01)l3-f 110)13)® 11)2 
:(|01)l3 + 110)13)® 10)2 



and the results for the concurrences are consistent with 
this form. C12 and C23 are vanishing while C13 increases 
to one. The sum of the squared concurrences is needed 
for calculating Q 



E^: 



(7/2 - (A - 2)77)^ 



(25) 



To calculate the tangle T123 one needs the i-concurrences 



ICi - IC3 



IC.^ 



1 

2^2 



A 



(26) 

(27) 



From these equations and the results for the concurrence 
it follows with l|12|) that ri23 = 0. Therefore the total 
global entanglement Q is given mainly by the squares of 
the concurrences as it follows from equation H16|) , and its 
parameter dependence is shown in fig. 2. 



TABLE I: Eigensystem of the 3-qubit Hamiltonian 



£1 = AJ 
E2 = AJ 



Ei = - 2- 



IV1) = I111) 
\i^2)^ 1 000) 
iV^3) = -^(|oii)-|iio)) 

i7/.4) = -^(|001)-|100)) 



E, = -|(r, + A-2) It^s) = ^ (1011) - ^1101) + 1110) ; 



2v^ 

£6 = -|(r, + A-2) \^.) = ^ 



E7 = 


-iiv- 


-A + 2) 


^8 = 

with 


-i(v- 


-A + 2) 



J100)-||010) + i001)j 



\^s) 



VvVx 
VvVx 



(1100) 



floio) 



1001)) 



rj = v'12 + A(A-4) and x = ^ + A - 2 



B. The states \tpr), l^s) 

The states \1lj7) and \ips) have the same results for en- 
tanglement measurement as well. The concurrences are 
calculated to 



C12 — C23 — — 



C 



13 



r?2-H(A-2)r?' 



(28) 



In the limit A — > 00 the concurrences are vanishing, the 
states have product form. For the special case A = 1, 
I ■07) and 108) have the form of the W-state 

107)- ^(1011) + |101) + |110)) = |I^) 
108) = -^(1001) + 1010) + |100)) = |iy) 



iK-J- 




FIG. 2: Comparison of global entanglement Q and the results 
of the Fs optimization for the states \%l}z) and \'4>7)\ addition- 
ally the optimization of {F-j,)'^ + {F-^)'^ for the state \tp7) is 
plotted. 



and the three concurrences are identical. The sum of the 
squared concurrences is needed for further calculations 



>;c^- 


«(? 




2 







(29) 


(7?2 + 


(A- 


- 2)77)- 


i-concurrf 


^nces 

-12 = 












-23 = ICz 


\/2 V 


2 - 

+ 


A 


(30) 


IC2 


-13 = 


2x/2 








(31) 



/Ci_ 



follows with 1)12(1 that ri23 = 0. Again there is no genuine 
3-qubit entanglement as the r— measure indicates. In 
addition the Q— measure is calculated (shown in fig. [SJ. 
Again, the global entanglement Q is only a function of the 
sum of the squared concurrences, which follows directly 
from equation H16|l . 

In the next section we compare these results with the 
optimized inequalities. 



C. Comparison with Bell optimization 

In fig. [3 besides the global entanglement Q we plot 
the results of the Bell optimization as a function of the 
anisotropy parameter A. The left y-axis shows Q, the 
right one the results for the squared inequalities. For 
both pairs of states \ip5), IV'e) and I'ipj), IV's) the inequal- 
ity ^ is violated. This is a sufficient condition for en- 
tanglement. With the squared inequalities @ and © we 
can distinguish 2-qubit and 3-qubit entanglement. The 
states 1^5) and lipg) are for A > 1.03 3-qubit entangled. 
The states JV't) and \^s) show 3-qubit entanglement in 
the range < A < 2.97. The points that mark the tran- 
sition between 2-qubit and 3-qubit entanglement were 
determined by {Fs)^ + {F^)l = 8.01. 
These results are in accordance with the 3-qubit classifi- 
cation after Diir et al. 113. There is the so-called W-class 
of states which are 3-qubit entangled and the tangle is 0. 
But if we compare the course of the optimized F3 with 



the course of the global entanglement Q we find - up 
to an scaling factor - exact analogy. This indicates that 
an apparent 3-qubit entanglement is due to to the fact 
that the sum of the squared concurrences is larger than 
1, the global entanglement Q larger than 2/3, cp. H16|l. 
In the discussion of these parameter dependent states an 
interesting result is that although both states have no 
finite tangle they differ in the aspect of the strength of 
the 2-qubit entanglement as measured by the sum of the 
squares of the concurrence. As soon as this sum is larger 
than 1 then there seems to be a kind of effective 3-qubit 
entanglement which is measured by the optimized Ft, op- 
erator. This means, that besides the pure 3-qubit tangle 
one has to consider the "strength" of the 2-qubit total 
concurrence, which might effectively describe some indi- 
rect 3-qubit entanglement. (But not a genuine one as 
measured by the tangle.) This means in our interpreta- 
tion that 3-qubit states with tangle equals to are only 
2-qubit entangled. 

With our results we conclude that also the W-state is 
only 2-qubit entangled, because the sum of the squarred 
concurrences is equal to 4/9, cp. ^fl]- 
In fig. 121 we plotted additionally the optimization of the 
squared inequalities for the state \ip7). For A > 4 the 
optimization yields max{(i^3)^ -|- (F^)^} = 8.00 indicatig 
2-qubit entanglement. As one can see, this method yields 
the same sufficient conditions due to entanglement clas- 
sification, but less information due to entanglement mea- 
surement. 

In the following we will discuss the superposition of two 
states in order to create a state with a finite tangle and 
to find at the same time the optimized F3 operator. 



D. Superposition of \tp7) and [tps) 

In this part we will discuss the superposition of the 
degenerated states \^pi) and j-^g): 



1 



V2 



(IV'r) + IV's)) 




FIG. 3: Measurement of the superposition of It/jt) and IV's); 
left y-axis: global entanglement, sum of squared concurrences 
and tangle; right y-axis: optimized F3 operator. 



In the limit A — > oo we get a highly entangled state 



VI. PURE 4-QUBIT HEISENBERG STATES 



V2 



(IV-t) + IV'8)) 



V2 



(|010) + |101)) (32) 



The concurrences, the i-concurrences and the tangle can 
be calculated exactly. The expressions for the concur- 
rence and the tangle are quite long and because of sim- 
plicity we will discuss them only graphically. 
The calculation of Q gives the following result: 



Q 



12 + hff + (A - 2)77 



(33) 



The Hamiltonian of our 4-qubit system can be written 
as 

H = -ricflai + ct"2(t^) + -^{02(^4), (34) 

with the product (Ti(fj — (jf(Jj + cf cj -I- (^iO'j- The cou- 
pling Js between spin 2 and 4 attaches the 4'^ spin to 
the 3 qubits interacting homogeneously. One can easily 
determine the eigenenergies and states of the system. We 
find out that two states are of special interest and call 
them |(i>i) and |$2)- They are energetically degenerated 
and belong to a spin triplet. The abbreviations we use 
in the following parts, are given in table II. 



In fig. I^lwe have plotted the sum of the squared concur- 
rences, the global entanglement Q and the tangle r as a 
function of A on the left y-axis. The right y-axis shows 
the results for the optimization of F^, . 
This superposition of two degenerated states has finite 
tangle and the sum of the squared concurrences is never 
larger than 1. From fig. Oit is interesting to note that at 
least above A = 3 there is a perfect agreement between 
the tangle and the optimized F3 operator. Below A = 3 
a tendency of the 3-tangle is reproduced. To sum it up it 
can be said that the F3 structure shows tangle like results 
and therefore measures genuine 3-qubit entanglement. 
Altogether one should also note that now the global en- 
tanglement measure Q sums up all different kinds of en- 
tanglement, the 2-qubit entanglement as measured by the 
sum of the squared concurrences and the 3-qubit entan- 
glement as measured by the tangle, cp. ((TC|l . 
It is remarkable in this respect that the F3 operator in 
its optimized form has an interesting structure. E.g. for 
large A where the state is mainly a GHZ type state, 
the operator mainly consists of linear combinations of a^ 
and (T^. There is no contribution from the ct^ compo- 
nents, but it is important to note that although the a^ 
contributions are small, they are important in the de- 
scription of the actual entanglement. This is seen by the 
fact that if one decreases A the entanglement decreases 
and this is seen by the fact that now the a^ contributions 
get much stronger although the a^ components still are 
negligible. Only for smaller A, namely in the region of 
A = 3 where the tangle goes to 0, one clearly sees that 
our optimized operator has now quite large contributions 
from a^. Looking into the calculation of the tangle one 
can conclude from this that a finite contribution in F3 
coming from the a^ operators may indicate a small gen- 
uine 3-qubit entanglement. 

This gives us sufficient confidence to discuss now 4-qubit 
systems where an explicit measure of 3- and 4-qubit en- 
tanglement is not known. And it turns out that the F4 
optimization will yield additional information. In order 
to compare with the 3-qubit results we use this time a 
special isotropic system (A = 1) and couple the 4**^ spin 
with a different coupling constant Jg. 



A. The state l$i) 

The state |<i>i) is of generalized W form and written as 

\^i) ^ aillllO) + 61IIOII) + cilOlll) - cilllOl). (35) 

Because of this structure it is clear from equation H18|) 
that the entanglement of this state is completely de- 
scribed by 2-qubit concurrences. These concurrences 
have been calculated in the following form: 



C12 — C'23 
C'la 

Ci4 — C34 

C24 



2V2 
1 



1 /^^8J(J-J,) 



-5 J + 2Js 



<52 



4J2 -J + 2Js 



1 



2\/2 V (52 

3 J + 2Js + 5 

AS 
iJ-Js){3 + 



9J-2Js 



6J + 2S 



(36) 

(37) 
(38) 
(39) 



Their dependences on the parameters J and Js are given 
in fig. 01i,b. With these results it is easy to sum the 
squares with the result given by 



j:c%'\{^ 



^2~ 



'J + 2Js 



(40) 



and the calculation of the global entanglement (cp. (|18|l 'l 
yields (sec fig. 0J:,d) 



1/ 12 J2 -J + 2Js 
= - 5- 



,52 



TABLE 11: Abbreviations for the 4-qubit model 

5 = V9 J2 - 4 JJ, + 4 J| 

1 



(41) 



0.1 = 
Cl = 



M2(3J+2Ja+a) 
2.y/l8J2-SJJs+8j2 



M2 - 

61 = 


\/3+ 3 

,/a-J 


V^M2{3J+5) 




1, {-J+2J,+S)2 




\l ^ ' 2.J--S 



1.0 



5 0.6 
gO.4 

u 

0.2 




a 

S 0.8 

B * 

a 

I 0.4 
■o 0.2 

o 
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1 1 1 1 1 1 1 

- (a) 
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FIG. 4: Concurrences, global entanglement and F4 optimization for the state |<l?i); (a): Cij{J,Js = 2); (b): Cij{Js,J ~ 2); 
(c): left y-axis: Q, right y-axis: F4 optimization, as a function of J, Js = 2; (d): left y-axis: Q, right y-axis: F4 optimization, 
as a function of Jj,, J — 2. 



From these figures various conclusions can be drawn. 
First of all, there are special points in the parameter 
space where certain concurrences are 0, especially for 
J = Js = 2. Form this one can conclude that the qubit 
2 can be separated in the state which indeed is true. 
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^(|011)i34 + |101)l34 - 21110)134) ® 11)2 
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Furthermore, it can be seen that a total entanglement 
decreases as a function of J and increases monotonically 
as a function of Js which means that by coupling these 
3 qubits to a 4'^ one in this special state we can increase 
the total entanglement which can be helpful in cetain ex- 
perimental situations. But most interestingly when we 
compare these results with the Bell inequality result, we 
find that the optimized Mermin-Klyshko polynomial op- 
erator F4 is up to a scaling factor the same function as the 
total global entanglement. But differently from this fac- 
tor in the case of F^ one can extract the information that 
for Js ^ 2 and J < 1.94 as well as for J = 2 and Js > 2.06 
the sum of the quadratic concurrences lies above 1 which 
equals to the fact that we have (Fi)^ + {FD^ > 8 which 
indicates an effective 3-qubit entanglement as described 
above due to the large sum of the squared concurrences. 
Another interesting aspect can be observed when the in- 
teraction constants are 0. When looking at the fig. ^ 
a,b all the 2-qubit concurrences are unequal to for 
Js — 2,J — and J = 2,Js — 0. At these values the 
entanglement is due to symmetry effects and not arising 
from interaction. 



B. The state |$2) 

Even more interesting are the entanglement character- 
istics fot the state |$2)- We will apply our reasoning also 
to this state although one can only give partially quanti- 
tative answers. First of all we note that the general form 
of the state given as 
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-1 1001) 



02 



|1100) 



-^10101) + ^11010) (42) 

is quite complicated but it reduces to a GHZ state in 
the two limits Js = 2 and J ^ as well as J = 2 and 

J, — !■ 00: 



1*2 
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- — (10101) -11010) 



(43) 



From this we can conclude that there must be besides 

the 2-qubit concurrences an additional 3- and/or 4-qubit 

entanglement. 

For the concurrences one finds: 
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(^12,^14,(723 and C34 are greater or equal than for 
Js — 2,J > resp. J = 2, Js > 0. The exact ana- 
lytic representation of C13 and C24 is only possible in 
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FIG. 5: Comparison of the sum of the squared concurrences and the F4 optimization for the state |$2). 



TABLE III: Parameters out of the F4 optimization for the state |$2); 
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the parameter ranges J — 2,Js < 2 and Jg — 2, J > 2: 

C'la — C2i — 

maxjo, ^(v/(52 + 5{J - 2Js) - U^- 
< \/2fi 



V52-J(J-2J,)-4J2)} (45) 



For J = 2, Js > 1 and Js = 2, J < 4, C13 and C24 are 
equal to 0. 

In fig. El we have plotted the sum of the squares of the 
concurrences. We can see that the maximum is at J^ — 
J = 2. It drops for very large Jg to while for J to 
infinity it levels to a finite value. Furthermore, it is found 
that the total global entanglement is constantly 1 



0(1*2)) 



1, 



(46) 



which means that there is no differentiation between the 
different qubit entanglements in this measure. Again it 
is very remarkable that the Mermin-Klyshko optimized 



operator is describing the additional entanglement (be- 
sides the 2-qubit concurrences) and follows parallel to the 
curve 1 — J2 ^ij ■ We therefore conclude in analogy to the 
3-qubit case that F4 measures the true 3- and 4-qubit en- 
tanglements for this state. This is the most interesting 
result of our paper, since here is a quantitative measure 
of n-qubit entanglement {n — 3,4) for a 4-qubit state, 
although we cannot discriminate between 3- and 4-qubit 
entanglement. 

Again if we look for the structure of F4 we find that all 
polynomial contributions for the limits J, = 2, J = and 
J ~ 2, Js —^ 00 come from products of cr^ and a^ for the 
different qubits and that the weight of this contribution 
changes with the strength of the additional entanglement 
(cp. Table Im)) . 

In a future paper the polynomial structure is investigated 
in detail. 



VII. CONCLUSIONS AND DISCUSSIONS 

It is shown in this paper that the optimized Mermin- 
Klyshko operators can be used very effectively to describe 
the degree of entanglement in different clusters of Heisen- 
berg spins. In those cases where there is in the 3-qubit 
system besides the concurrences no additional entangle- 
ment (i. e. the tangle r = 0) the optimized F3 op- 
erator perfectly describes the 2-qubit entanglement of 
the system as a function of the anisotropic parameter 
in the Heisenberg cluster and it is more or less identical 
to the global entanglement measure resp. the sum of the 
squared concurrences, cp. Ijltill . In those cases where in 
addition to the 2-qubit concurrences there is a finite tan- 
gle r, we find that this additional 3-qubit entanglement 
measured by r is nearly perfect described by the opti- 
mized F3, as shown in fig. O 

We therefore test this result in a 4-qubit system and 
again we find two different cases. We discuss a special 
state l^i) where the sum of the 2-qubit concurrences is 
mainly proportional to the global entanglement measure 
and from this we find that the optimized F^ operator 
follows this Q value. In the second eigenstate for the 



system |<I>2) where the global entanglement is just equal 
to 1, independent of the parameters, we expect besides 
the 2-qubit concurrences an additional entanglement and 
this seems to be perfectly the case, especially when the 
expression 1 — X) ^ij i^ compared in its parameter de- 
pendence to the optimized -F4, shown in fig. [3 Our 
special interest is here that even at the minima of these 
functions, at the point J = Jg = 2, there is a small but fi- 
nite higher entanglement which of course at the moment 
could not be interpreted as 3- or otherwise 4-qubit en- 
tanglement. 

It should be noted that the optimization procedure for 
the F]s[ operators heavily depends on the starting values 
and therefore a procedure has to be used where a random 
choice for the starting values has to be done. Another 
remark is, that this optimization yields much more than 
one minimum or maximum and therefore one should be 
careful with the interpretation of these parameters. But 
at least for the GHZ-state with 4-qubits we have shown 
that these operators contain besides the usually used a^ 
operators additional a^ contributions \i^ . 
Further work is in preparation where a more extensive 
study of these F]\j operators will be presented. 
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